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We present a qualitative (and quantitative, at the level of 
estimates) analysis of the ordering kinetics in a strongly non- 
equilibrium state of a weakly interacting Bose gas, trapped 
with an external potential. At certain conditions, the or- 
dering process is predicted to be even more rich than in the 
homogeneous case. Like in the homogeneous case, the most 
characteristic feature of the full-scale non-equilibrium process 
is the formation of superfluid turbulence. 

PACS numbers: 03.75.Fi, 32.80.Pj, 67.40.-w 

Kinetics of Bose-Einstein condensation (BEC) in a 
weakly interacting Bose gas is one of the most funda- 
mental problems of non-equilibrium statistical mechan- 
ics. The exciting progress in the experiment on BEC in 
ultracold gases initiated by the pioneer works |ij opens 
up an opportunity of laboratory study of BEC kinetics. 
Though an important step in this direction has been al- 
ready made (2l , most non-trivial regimes of ordering ki- 
netics have not been achieved yet. 

In Refs. &fl] (see also Ref. H ) a scenario of BEC in a 
strongly non-equilibrium self-evolving homogeneous gas 
was developed. The scenario involves a number of qual- 
itatively different stages. The evolution starts with an 
explosion-like wave in the wavenumber space, propagat- 
ing from higher energies towards the lower ones, and re- 
sulting in a formation of a specific power-law distribution 
of particles over the wavenumbers [p| . Immediately after 
its formation, this distribution starts to relax, the relax- 
ation having the form of a back wave in the wavenumber 
space |j. Simultaneously, in the low-energy region the 
so-called coherent regime sets in that leads to the forma- 
tion of the quasicondensate state Q . Basically, the quasi- 
condensate state corresponds to what is known in the the- 
ory of superfluidity as the state of superfluid turbulence. 
It can be viewed as a condensate containing a tangle of 
vortex lines (plus a sharply non-equilibrium distribution 
of long- wave phonons) B. The formation of the qua- 
sicondensate occurs very rapidly (characteristic time is 
much smaller than the time of the wave formation) p,EJ . 
In contrast to it, the final stage of long-range ordering, 
associated with relaxing superfluid turbulence and long- 
wave phonons, takes a macroscopically large (the size of 
the system enters the answer) time S . 

From a common wisdom one may expect that at 
certain conditions strongly non-equilibrium BEC in a 
trapped gas should reproduce main qualitative features of 
the homogeneous scenario. It is not surprising thus that 
the explosion-like initial kinetic stage [pj is quite relevant 
to the case of the trapped gas fj\M and is consistent with 



the experimental observations of Ref. M. Nevertheless, 
full-scale experimental realization of the scenario y|§] in 
a trapped gas is not guaranteed yet by just considerable 
deviation from the equilibrium in the initial state. The 
point is that the situation with the BEC kinetics in ultra- 
cold gases is determined by a competition between large 
number of particles and very small interaction. The lat- 
ter circumstance leads to a rather large coherent length 
associated with the ordering process, which for a trapped 
gas can turn out to be just larger than the size of the con- 
densate, in which case the character of the condensate 
formation will be drastically changed: Instead of qua- 
sicondensate state there will take place a build-up of a 
genuine condensate with the initial form pre-determincd 
by the external potential. The theoretical analysis of the 
BEC kinetics in a trapped gas has been dealing mostly 
with this physical picture |7|q|, ignoring thus the non- 
trivial regime associated with the superfluid turbulence 
formation. 

The main goal of the present letter is to consider the 
superfluid-turbulcncc regime of BEC kinetics in an ex- 
ternal potential. We will demonstrate that the necessary 
and sufficient condition for this regime to occur corre- 
sponds to violation of the Knudsen condition at a certain 
stage of self-evolution within a certain region of coordi- 
nate space. While the external potential together with 
the initial number and energy of particles determines the 
initial paranreters of the quasicondensate: the size i?*, 
the number of quasicondensate particles Nq, and typical 
initial separation between vortex lines in the vortex tan- 
gle, rg , the very process of the quasicondensate formation 
goes in a quasi-homogeneous anti-Krmdsen regime when 
the external potential becomes irrelevant and, in prin- 
ciple, may be even turned off without any considerable 
effect on the process of the quasicondensate formation. 

In the experiments with trapped ultracold gases nor- 
mally there takes place the Knudsen regime, when the 
free path length of a particle with the energy e, lf rC c{s, n e ) 
[which depends also on occupation numbers n e because 
of the stimulated scattering], is much larger than the 
typical radius of the particle's trajectory, R e . For def- 
initeness we consider a parabolic trap with all the three 
frequencies of the same order c^o- So that R £ ~ loov(e) 
[v(e) is the typical velocity corresponding to the energy 
e], and the condition lf lee (e,n £ ) 3> R e is equivalent to 
T ou(e)»e)wo > 1, where (Ti = 1) 



T co n(s,n e ) ~ l/m(asn e y 



(1) 



is the collisional time. Here a is the scattering length and 
m is the particle mass. 



Knudsen regime is a convenient starting point for an- 
alyzing kinetics in a potential. Almost in all qualitative 
aspects it corresponds to an isotropic homogeneous case, 
since the distribution of particles depends only on the two 
variables, e and t (ergodic approximation). The scaling 
of the collision time remains the same, so that the main 
quantitative difference comes from the difference in the 
density of states. 

During some initial period of evolution the correlations 
between different single-particle modes are vanishingly 
small. Such a regime admits a description in terms of 
kinetic equation. This stage is thus referred to as kinetic 
stage. If the deviation from equilibrium is strong enough, 
the evolution at the kinetic stage results in the formation 
of a self-similar wave in the wavenumber space propa- 
gating in an explosion-like fashion from the high-energy 
region (where the particles are initially distributed) to- 
wards lower energy scales. Corresponding solution of the 
kinetic equation has the form p] 



n £ (t) = Ae^ a (t)f(s/e (t)), t<t. 



(2) 
(3) 



Here A and B are dimensional constants depending on 
the initial condition and related to each other by the 
formula B = const(ma 2 A 2 ) 1 / 2 ^ Q_1 ^. The dimensionless 
function / is defined up to an obvious scaling freedom. 
The explosion character of the evolution guarantees that 
the wave reaches the point e = at some finite time mo- 
ment t = t* [t = corresponds to the beginning of evolu- 
tion] , the value of i* being on the order of the collisional 
time at t = 0. Physically, the explosion-like evolution 
is supported by the stimulation of the collision rate at 
the head of the wave, e ~ £o, by ever growing occupa- 
tion numbers. A remarkable feature about the solution 
(0)"@ i s it s universality: The coefficient A is the only 
quantity related to a particular cooling regime or a par- 
ticular form of the initial condition. Having in mind the 
most characteristic statement of the problem when the 
self-evolution starts with an essentially non-equilibrium 
state of N particles in a parabolic trap of typical fre- 
quency ujq, with a typical occupation numbers 7 > 1, we 
readily estimate A by extrapolating the solution (0)-(||) 
to t = 0: 
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Generally speaking, the index a in (g)-(^) cannot be es- 
tablished from the scaling properties of the collision term 
of the kinetic equation, being related thus to the particu- 
lar form of the latter M. It is possible, however, to specify 
lower and upper limits for a following from the consis- 
tency of (0)-(0) with the requirement that these formulae 
describe an explosion-like singularization of distribution 
(rather than infinite-time shrinking). To this end we note 
that from the scale invariance of the collisional term it 



follows that f(x) should behave at x ^> 1 like some power 
of x. Then, from the fact that at t = £* the occupation 
numbers are finite at e > one concludes that 



m 



at x 



(5) 



The requirement that the particle distribution does not 
shrink as a whole implies that the number-of-particles 
integral for the distribution (||), (||) is divergent at e — > 
00. The number of particles N = J w(e)n e de , where the 
density of states w(e) is given by 

( \ _ I Fm 3 / 2 e 1 / 2 /v / 27r 2 (homogeneous) , ,.,. 

1 2e 2 /w x uj v w z (parabolic potential) . 

Here V is the volume of a homogeneous system and ui x , 
uiy, and u z are the frequencies of the 3D harmonic oscil- 
lator. This immediately yields a < 3/2 for the homoge- 
neous case and a < 3 for the case of parabolic potential. 
The condition a > 1 is necessary for £o(*) to approach 
zero at t = t*. So we have 1 < a < 3/2 in the homoge- 
neous case and I < a < 3 in the parabolic potential. The 
most accurate up-to-date numeric analysis of a in the ho- 
mogeneous case was performed in Ref. H with the result 
a « 1.24. In the case of parabolic potential a rs 1.6 |[Lp|. 
Formally, at t = t* the front of the wave £o(i) reaches 
the energy e = 0. However, somewhat earlier the solu- 
tion (|2])-([|) becomes inapplicable in the region of very 
small energies. In the homogeneous case this occurs due 
to the change of the evolution regime (in corresponding 
low-energy region) from weak turbulence to strong tur- 
bulence (the so-called coherent regime) P,B- In the ex- 
ternal potential, however, the scenario is modified. As we 
demonstrate below, what changes the solution (0)-(0) at 
a certain stage of evolution of a trapped gas is either (i) 
discreteness of the levels, or (ii) violation of the Knud- 
sen regime in a certain region of the coordinate space. 
Apart from a cross-over region, the two circumstances 
are mutually exclusive: If discreteness of the levels be- 
comes relevant, the Knudsen regime will not be violated, 
and if the Knudsen regime is violated, the discreteness 
of levels will not become relevant. Hence, the change of 
the evolution in a trapped gas from the self-similar wave 
(g)-(||) is governed by the following two parameters: 

£ = Tcoii^o (Knudsen parameter) , (7) 



d = T co \\/w{e) (discreteness parameter) . (8) 

Knudsen regime corresponds to £ ^> I; the discreteness 
of the levels is irrelevant when d« 1. Note that r co ii 
essentially depends on s and n £ , which ultimately means 
that both £ and d are the functions of energy and time. 
First, let us consider the case (i), where the condition 
d ~ 1 occurs when the front of the wave £o(£) reaches 
some energy e ( i- As it follows from ([!])- (||), (||), the pa- 
rameter d at the front of the wave can be estimated as 



3 2(e*-2) 
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(9) 



From Eq. (||) we see that d So decreases with decreasing £o 
only if a < 2, which means that the case (i) is possible 
only under this condition. [In this respect the result a sw 
1.6 MM is of qualitative importance!] Substituting A 



from 



into (P|) and setting d eo ~ 1, we find 



Sd ~ w p 



-1/2(2- 



p = uj Q ma' 



(n° 



/3 7 1" 



■a/3 



(10) 



(11) 



Here p is the main dimensionless parameter which we 
will encounter several times. By definition, in the case 
(i) the value of £d Eq. (|l0|) should be greater than luo, 
which implies p« 1. Below we will see that p -C 1 is 
not only a necessary, but also a sufficient condition for 
the case (i) to occur: Under this condition the case (ii) 
is impossible. 

To describe what happens after the front £o(t) reaches 
the energy scale e^, we need to make an assumption that 
the discreteness of the levels leads to a slowing down of 
the evolution cascade at e < £d- Such an assumption ]1C| ] 
seems to be quite natural in view of the fact that discrete 
harmonics of different frequencies are out of resonance, 
which should suppress the cascade in the wavenumbcr 
space. If this is the case, the next stage of evolution at 
the energies £ < £d will be a fast equilibration [with the 
characteristic time r co ii(s = £d,n e = A/e°[); harmonics 
with energies e <C Ed equilibrate due to the interaction 
with harmonics e ~ Ed, rather than with each other]. 
Further evolution should lead to the destruction of the 
non-equilibrium distribution n £ = A/e a at e > £^, the 
process having the form of a self-similar back cascade in 
the wavenumber space (back wave) |3[] : 

n e (t) = AEo a (t)f(E/s (t)) , e>0, t>U (12) 

[f(x) — ► f(x) at x — » oo], where £o(t) is given by the 
same formula (|j). Clearly, Eq. (^2|) is applicable only 
if t — i* greater than the above-mentioned time of equi- 
libration of the harmonics e < Ed- As follows from an 
estimate of corresponding terms of collision integral || , 
the back wave should create a quasi-equilibrium distri- 
bution at e <C £a{t), that is f(x) oc 1/x at x -C 1, with 
infinite at t = £* and ever decreasing afterwards temper- 
ature oc £Q a (t). 

Since after the formation of the back wave all the har- 
monics with e <C £o{t), including the lowest one (!), are in 
a quasi-equilibrium, the process of long-range ordering in 
the case (i) is somewhat trivial: The state of the harmon- 
ics £ <C £o(t) is determined by the instant temperature. 
In particular, the number of condensate particles No(t) 
follows from Eq. (O) by conservation of the total number 
of particles: 



N (t) 



2.4 
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■£^ a (t) I dxx 2 \x 



/»]• (13) 



This universal law of the condensate build-up should take 
place until the back wave cascade reaches the initial en- 
ergy scale and the further growth of Nq becomes sensitive 
to the details of preparation of the initial state. 

We now turn to the most non-trivial case (ii), where 
the condition £ ~ 1 occurs when the front of the wave 
£a(t) reaches some energy eg. From (0)-(||) and (g) we 
have the following estimate for the parameter £ at the 
front of the wave: 

£ £0 ~ we^-V /ma 2 A 2 . (14) 

Substituting A from (||) into (jlj) and setting £ eo ~ 1, we 
find 



£e 



uj p 



l/2(a-l) 



(15) 



which, in particular, means that the case (ii) can only 
occur if p > 1. From Eqs. ( fl0| ) and ([F|) we conclude 
that, apart a cross-over region (corresponding to p ~ 1), 
the following alternative takes place. Either discreteness 
of the low- lying levels becomes relevant (p <C 1), or there 
takes place a break-down of Knudsen regime (p ^> 1). 
Geometric size i?* of the area where the break-down of 
Knudsen regime takes place is obviously related to e^ 
by muj^R 2 ~ e^, which gives i?* ~ ZoP 1 ^ 4 ^ -1 ^ where 
Iq = 1/y/mujo is the typical size of the single-particle 
groundstate wavefunction in the trap. 

A natural question then is: What happens at distances 
r < R* from the center of the trap after the Knudsen 
regime is broken down? A self-consistent scenario of fu- 
ture evolution can be constructed by noticing that fur- 
ther increase of the collisional rate due to the stimulated 
scattering naturally leads to the anfi-Knudsen regime, 
when the free-path-length of the low-energy particles be- 
comes much smaller than i?». This immediately leads 
to a quasi-homogeneous picture of evolution: the region 
r < R* can be formally thought of as a number of inde- 
pendent homogeneous sells of the size much larger than 
the free-path- length, but much smaller than i?*. Evolu- 
tion in each sell follows homogeneous scenario, ultimately 
leading to the formation of the quasicondensate. Indeed, 
in the anti-Knudsen regime the evolution during the time 
period on the order of collision time is insensitive to the 
external potential (the criteria for the Knudsen regime 
and for the sensitivity to the potential within the colli- 
sion time coincide). But this time is enough to form the 
wave (0)-(0) in the wavenumber space (with the expo- 
nent a now corresponding to the homogeneous case) and 
then to form quasicondensate. During the evolution in 
the wavenumber space, the free-path length of the parti- 
cles with £ ~ £o(i) is getting progressively smaller, which 
renders the proposed scenario self-consistent. 

A minor deviation from the pure homogeneous picture 
is that now the moment i* depends on the distance from 



the center of the potential, so that the coherent regime 
first should start at r = (the point of maximal initial 
density) and then gradually occupy all the anti-Knudsen 
region up to r ~ R*. Thus, the process of quasiconden- 
sation takes on a form of wave propagating in the coordi- 
nate space from r = towards r ~ R*. When this wave 
reaches r ~ i?* , the process of formation of the quasicon- 
densate correlation properties is finished, and within the 
region r < i?* we have a state of superfluid turbulence. 

Note that the fact that the quasicondensate forma- 
tion in the anti-Knudsen regime goes independently of 
the external potential should lead to certain qualitative 
effects. For example, after the anti-Knudsen regime is 
well formed, the trapping potential can, in principle, be 
turned off without any dramatic effect on the (short- 
range) ordering process. Clearly experimental realiza- 
tion of such a regime would be interesting and instructive 
from the fundamental point of view. Another effect, com- 
ing from the same reason, is that the quasicondensate, 
just upon its formation, is out of equilibrium with re- 
spect to a global coherent motion in the potential, which 
means that a breathing mode for the quasicondensate 
motion should ultimately be excited in a trapped gas. 

As far as the evolution in the anti-Knudsen region is 
of quasi-homogeneous character, all the quantities char- 
acterizing this process can be estimated from homoge- 
neous relations of Refs. 0-0]. As an initial condition to 
the corresponding homogeneous case one should take an 
extrapolation of the inhomogeneous solution (0)-(0) to 
the region of breaking down the Knudsen regime. Hence, 
the typical initial homogeneous single-particle energy and 
occupation number are nothing else than s^ and A/ef, 
respectively. Correspondingly, the number of particles 
involved in the quasi-homogeneous BEC process can be 
estimated as 



R*/r ~p v , v 



4a n 



4(a-l)(2a -l) 



(17) 



N ^ y-a/3jya/3p(3-c0/2(a-l) 



(16) 



An important step of the qusi-homogeneous regime is as- 
sociated with coming the homogeneous back wave (§)-@ 
to the initial energy scale (which in our case corresponds 
to £<;). At this stage the number of qussicondensate par- 
ticles Ao is on the order of the total number of particles. 
Eq. ( |l6| ) thus yields an estimate of the number of quasi- 
condensate particles which will be created in the course 
of anti-Knudsen self-evolution in the system, even if the 
trapping potential is turned off upon the formation of the 
anti-Knudsen regime. 

The state of superfluid turbulence, which arises to- 
gether with the formation of the quasi-condensate local 
correlation properties, is characterized by the typical sep- 
aration between the vortex lines. This separation ro is 
minimal just upon the quasicondensate formation, when 
ro is on the order of the quasicondensate healing length. 
From the estimate of the quasicondensate density upon 
its formation 11 we find ro and come to the following 
relation: 



where ao rs 1.24 is the exponent a for the homogeneous 
case. From (|l7| ) we see once again that the formation of 
the superfluid turbulence is possible only if p 3> 1 . 

In conclusion, we note that in the experiment of Ref. 
g the parameter p is, roughly speaking, of order unity. 
This means that the full-scale non-equilibrium BEC sce- 
nario cannot be seen with these initial conditions. The- 
oretically, the most direct way of achieving p 3> 1 is to 
radically increase ujq immediately after the formation of 
the non-equilibrium initial state, during the time which 
is much less than the collisional time. 
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In this connection we cannot but note that the un- 
dertaken in Ref. M attempt of a priori identifying a 
with the exponent of Zakharov's constant-particle-flux 
steady-state solution of the kinetic equation is inade- 
quate. (Which, in particular, was demonstrated by an 
accurate numeric analysis of Ref. |6J.) 
An accurate proof of this assumption, as well as a quanti- 
tative check of the other our qualitative predictions could, 
in principle, be done by numerically simulating time- 
dependent Gross-Pitaevskii equation (GPE) with corre- 
sponding turbulent initial conditions, the applicability of 
GPE being guaranteed by large occupation numbers [see, 
e.g., Yu. Kagan and B.V. Svistunov, Phys. Rev. Lett. 79, 
3331 (1997)]. 



